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Abstract — Due to the increasing demand of capacity in wireless 
cellular networks, the small cells such as pico and femto cells are 
becoming more popular to enjoy a spatial reuse gain, and thus cells 
with different sizes are expected to coexist in a complex manner. 
In such a heterogeneous environment, the role of interference 
management (IM) becomes of more importance, but technical 
challenges also increase, since the number of cell-edge users, 
suffering from severe interference from the neighboring cells, will 
naturally grow. In order to overcome low performance and/or high 
complexity of existing static and other dynamic IM algorithms, 
we propose a novel low-complex and fully distributed IM scheme, 
called REFIM (REFerence based Interference Management), in 
the downlink of heterogeneous multi-cell networks. We first 
formulate a general optimization problem that turns out to 
require intractable computation complexity for global optimality. 
To have a practical solution with low computational and signaling 
overhead, which is crucial for low-cost small-cell solutions, e.g., 
femto cells, in REFIM, we decompose it into per-BS (base station) 
problems based on the notion of reference user and reduce 
feedback overhead over backhauls both temporally and spatially. 
We evaluate REFIM through extensive simulations under various 
configurations, including the scenarios from a real deployment of 
BSs. We show that, compared to the schemes without IM, REFIM 
can yield more than 40% throughput improvement of cell-edge 
users while increasing the overall performance by 10~107%. This 
is equal to about 95% performance of the existing centralized IM 
algorithm (MC-IIWF) that is known to be near-optimal but hard 
to implement in practice due to prohibitive complexity. We also 
present that as long as interference is managed well, the spectrum 
sharing policy can outperform the best spectrum splitting policy 
where the number of subchannels is optimally divided between 
macro and femto cells. 

Index Terms — Interference management, heterogeneous wireless 
access networks, femto cells, reference user, power control, user 
scheduling, feedback reduction, distributed algorithm; 



I. Introduction 

Many researchers from networking and financial sectors 
forecast that by 2014, the total mobile data traffic throughout 
the world will grow exponentially and reach about 3.6 exabytes 
per month, 39 times increase from 2009 iQ], Pushed by this 
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explosive demand mainly from bandwidth-hungry multimedia 
and Internet-related services in broadband wireless cellular 
networks, communication engineers seek to maximally exploit 
the spectral resources in all available dimensions. 

Small cells such as pico and femto cells, seem to be one 
of the most viable and economic solutions |[3]. Of recent 
significant interest is the femto cell designed for usage in 
a home or an office and deployed by users, utilizing user's 
Internet connection, e.g., cable or DSL (Digital Subscriber 
Line) Internet services as a backhaul. Macro and pico cells 
are controlled by mobile network operators and use dedi- 
cated backhauls. Small cells are also considered as a way of 
incrementally increasing coverage and/or capacity inside the 
initial deployment of macro cells. In addition to the advantages 
of spectrum reuse efficiency, the network operators are also 
attracted by financial benefits because small cells can reduce 
both capital (e.g., hardware) and operating (e.g., electricity, site 
lease and backhaul) expenditures. 

In heterogeneous multi-cell networks with a mixture of 
macro and small cells, interference is a major obstacle that 
can impair the potential gain of small cells and its pattern is 
highly diverse [4J, e.g., interferences in macro-to-macro, femto- 
to-femto, and macro-to-femto. As the number of small cells 
increases, the number of users at cell edges suffering from low 
throughput due to severe interference also grows. In particular, 
femto base stations (BSs) are installed in an ad-hoc manner 
without being planned by users, not the network operators, 
which also increases the technical challenges of interference 
management (IM). 

To mitigate interference, a traditional frequency reuse or 
more enhanced schemes such as fractional frequency reuse 
(FFR) lis and its variations |l6l, Q can be utilized. All these 
schemes represent static IM algorithms for pure macro cell 
networks, where a specific reuse pattern is determined a priori 
by the network operator at offline. However, in reality, BSs are 
not uniformly deployed over the network. In particular, femto 
BSs are purely controlled by users, which implies that they may 
often be installed by users and even existing ones may be turned 
on/off dynamically. Under this situation, having the static reuse 
pattern is naturally inefficient. Recently, several dynamic IM 
algorithms have been proposed to address this problem [8 |- 
|14|. They can significantly improve the performance over the 
static schemes, but many of them suffer from prohibitively high 
complexity and message passing among neighboring BSs. 

In a single-carrier multi-cell network setup, several BS 
coordination schemes ISl- lfTol have been proposed under the 
assumption of binary power control, i.e., each BS transmits data 
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with its given maximum power or zero. There are several recent 
works in a multi-carrier multi-cell network JTT], ^T\- Venturino 
et al. [ 1 1 1 proposed several algorithms which perform multiple 
iteration loops in a slot for user scheduling and power allo- 
cation. Although they can achieve near-optimal performance, 
all of them are centralized algorithms which are too complex 
to be implemented in practice. Stolyar et al. ||T21 proposed 
algorithms that adjust BS powers much more slowly than per- 
slot user scheduling. This time-scale separation does simplify 
the problem solution and reduce the complexity, but may lead 
to non-negligible performance loss. To tackle the cross-tier in- 
terference between macro and femto cells, there also have been 
several approaches by Sprint, Ericson llT3l and Chandrasekhar 
et al. U?] that adjust the transmit powers of femto BSs based on 
the relative locations with respect to the macro BS. The power 
control also has been handled in ad-hoc networks. Chiang et al. 
ifTSi showed that in high SINR (signal to interference plus noise 
ratio) regime nonconvex power control optimization problems 
can be transformed into convex optimization problems through 
a geometric programming technique. In addition, there has been 
research on mitigation interference in a MIMO (multiple input 
multiple output) setting 1161 . 

The part of dynamic IM algorithms, in particular, the power 
control component, was studied in the context of multi-tone 
DSL networks (see ifTTl and the references therein for a nice 
survey). Indeed, the wired multi-tone DSL model with crosstalk 
can be interpreted as a special case of the wireless multi-carrier 
cellular model with inter-cell interference when (i) only one 
user exists per cell and (ii) wireless channels are stationary 
and (iii) distributed operations among BSs are not crucial. In 
fact, the power control of our proposed solution is motivated by 
ASB (Autonomous Spectmm Balancing) HT), GSjj in the DSL 
network that uses the idea of reference line. However, in the 
IM over multi-cell wireless networks, much more challenging 
issues still remain for practical implementation, e.g., joint 
operation with multi-user scheduling in each cell, dynamic 
selection of reference users over time-varying channels, and 
small message passing among neighboring cells. 

The fundamental challenges in the dynamic IM are that 
(i) BS power control problem itself (even if scheduled users 
in each BS are fixed) is formulated by a highly nonconvex 
optimization ifTTI . lfT2ll . ifTSI, it is tightly coupled with multi- 
user scheduling, and heavy message passing is usually 
required to coordinate BS powers. In this paper, we aim at 
developing an IM scheme consisting of joint power alloca- 
tion and user scheduling which is practical in terms of low 
complexity and small message passing, but yet the scheme 
achieves near-optimal performance. Low complexity and small 
message passing is particularly essential for low-cost solutions 
such as femto BSs because they are typically made of cheap 
devices for price competitiveness and connected to the low- 
speed residential cable or DSL Internet connections, not to the 
high-speed dedicated backhauls. 

Another important issue in heterogeneous networks is the 
way of sharing spectrum between macro and femto cells. If 
network operators adopt a spectrum splitting policy where 
macro and femto cells orthogonally use the resource for conve- 
nience of implementation, they can be free from the macro-to- 



femto interference. However, such a splitting policy needs to 
determine the ratio of optimal splitting that varies depending on 
the configuration of cells, e.g., density and position, resulting in 
spectrum inefficiency. Alternatively, a spectrum sharing policy 
between macro and femto cells can be adopted to maximally 
reuse the resource. However, the macro-to-femto interference 
may harm the system performance unless appropriate IM al- 
gorithms are employed. It will be interesting, especially for 
the network operators, to answer which of these two polices is 
better under what circumstances. 

Motivated by the above, we propose a novel IM algorithm, 
called REFIM (REFerence based Interference Management), 
whose core features are summarized as follows. 

1) With the notion of reference user, each BS can approximate 
the interference impact of all other cells with a single user 
to which the BS generates the most significant interference. 
This abstraction substantially simplifies the problem, result- 
ing in the power control algorithm with low computational 
and signaling overhead. 

2) Due to the nonconvexity of the power control problem, 
different initial power settings may lead to different solu- 
tions. We empirically show that running the power control 
algorithm with the powers used at the previous slot as the 
initial powers can have effect of removing multiple loops 
without much performance degradation. 

3) In the original power control based on the reference user, it 
requires to feedback per-user information in a cell to the 
neighboring BSs at every slot. We reduce such a heavy 
message passing overhead over backhauls both temporally 
and spatially. 

4) REFIM has a nice feature of incremental deployability 
that partial deployment in some specific regions, probably 
starting from the regions that experience small capacity due 
to severe interference, sufficiently increases the capacity in 
those regions, but not affecting other regions. This is a 
desirable property for network operators who cannot afford 
to upgrade the IM module in all BSs. 

5) We also demonstrate that, as long as an appropriate IM 
such as REFIM is adopted, the performance of the spectrum 
sharing policy is much better than the best performance of 
the spectrum splitting policy that the number of subchannels 
is optimally divided between macro and femto cells. 

The remainder of this paper is organized as follows. In 
Section ini we formally describe our system model and general 
problem. In Section InD we propose a reference user based 
power allocation and user scheduling with feedback reduction 
ideas. In Section IIVI we analyze the computational and sig- 
naling complexity. In Section [V] extensive simulations under 
various configurations demonstrate the performance of the 
proposed algorithm compared to previous algorithms. Finally, 
we conclude the paper in Section [VT] 

II. System Model and Problem Definition 

A. Network and Traffic Model 

We consider a wireless cellular network consisting of mul- 
tiple heterogeneous BSs, where Mmacro, J^pico and M femto 
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are the set of macro, pico and femto BSs, respectively. Denote 

by /C = {1,...,K} and TV = Mmacro U Npico U Nfevato = 

{1,...,A^} the set of users and BSs, respectively. BSs and 
users are equipped with one transmit and one receive antenna, 
respectively. Each user is assumed to be connected to a single 
BS. Denote by /C„ the (nonempty) set of users associated with 
the BS n, i.e., /C = /Ci U • • • U Kn and /C„ ("1 /C™ = 0, for 
n ^ m. A full buffer traffic model with infinite data packets 
in the queue for each user at its associated BS is used to 
consider best-effort traffic. Macro and pico BSs can in general 
exchange information very fast with each other because there 
are connections between them via high-speed wired dedicate 
backhauls directly or through a base station controller. On the 
other hand, femto BSs can be assumed to exchange information 
slowly through user's Internet connections as backhauls. 

B. Resource and Allocation Model 

We consider a system where a subchannel is a group of 
subcarriers as the basic unit of resource allocation. Assume 
that there are S number of subchannels and all BSs can use all 
the subchannels for data transmission, i.e., universal frequency 
reuse. Denote by S = {1, . . . , 5} the set of subchannels. We 
focus on the downlink transmissions in the time-slotted system. 
At each slot, each BS needs to determine (i) which user is 
scheduled on each subchannel and (ii) how much power is 
allocated for each scheduled user on each subchannel. 

User scheduling constraint: In regard to (i), denote by 
I^[t) = [I^'^{t) : fc e /C, n G A/] the user scheduling indicator 
vector, i.e., I^''^{t) ~ 1 when BS n schedules its associated 
user k on subchannel s at slot t, and otherwise. Furthermore, 
we denote the user scheduled by BS n on subchannel s at slot 
t by k{n,s,t). Reflecting that at most only one user can be 
selected in each subchannel for each BS, we should have: 

J2 Is'"{t) < 1, Vn e A/", s e S. (1) 

Power constraint: In regard to (ii), denote the transmit power 
of BS n on subchannel s at slot t by Pg{t). The vector 
containing transmit power of all BSs on subchannel s is 
Ps{t) = [pl{t), ■ ■ -p^ ■ In parallel, the vector containing 
transmit powers of all subchannels for BS n is = 
• • ■p's{i)V ■ Each BS is assumed to have the total power 
budget and spectral mask constraints: 

^p^(i) < P"''"'^^ VneAA, (2) 

< p;'''"'^^^ Vn e AA, s e S. (3) 

In practice, a typical transmit power of macro BSs is around 
43dBm, which is 20~30dBm higher than that of small BSs. For 
notational simplicity, the time-slot index it) is dropped unless 
confusion arises. 

C. Link Model 

In this paper, we do not consider advanced multiuser de- 
tection or interference cancellation, and hence the interference 
from other BSs is treated as noise. We focus on the spectrum 



level coordinatiorQ, i.e., finding multi-channel power allocation 
of each BS in order to improve system performance by mitigat- 
ing the interference. For a given power vector Ps, the received 
SINR for user k from BS n on subchannel s can be written as: 

7^'"(P.) - = 7, (4) 

where p" and " representing the nonnegative transmit power 
of BS n on subchannel s and the channel gain between BS n 
and user k on subchannel s during a slot, respectively; u^' is 
the noise power. The channel gain is time-varying and takes 
into account the path loss, log-normal shadowing, fast fading, 
etc. 

Following the Shannon's formula, the achievable data rate 
[in bps] for user k on subchannel s is given by: 

r^"(p,) = |log2(^l + i7''"(p.)), (5) 

where B denotes the system bandwidth; F denotes the SINR 
gap to capacity which is typically a function of the desired 
bit error ratio (BER), the coding gain and noise margin, e.g., 
F = - iii(oBER) M-QAM (quadrature amplitude modulation) 
1191 . Note that r^'^{ps) is the potential data rate when the user 
k is scheduled for service by BS n on subchannel s and its 
actual data rate becomes zero when another user is scheduled, 
i.e., rg'"(ps, Is) = Is'" ■ Ts'^iPs)- We assume that F = 1 and 
drop B/S mainly for simplicity, but our results can be readily 
extended to other values of F and B/S. 

D. General Problem Statement 

Our objective is to develop a slot-by-slot joint power 
allocation and user scheduling algorithm that determines 
{p{t))Z^ and {m)Z^, where p{t) = {ps{t),s e S) 
and I{t) = {Is{t),s £ S). The long-term achieved 
throughput vector R = {Rk : k G K,), where = 
^i^jJ2l=iJ2ses'^s''''iPsi'^)'^s{T)), is the solution of the 
following optimization problem: 

(Long-term P) : max ^ Uk{Rk) (6) 

keK 

subject to KeTZ, (7) 

where Uk{-) is a concave, strictly increasing, and continuously 
differentiable utility function for user fc; 7e C is the 
set of all achievable rate vectors over long-term, referred to 
as throughput region. Note that this optimization problem is 
challenging to solve, since we aim at devising an instantaneous, 
distributed algorithm even though the constraint set TZ is neither 
available to the users nor the BSs. 

With the help of the stochastic gradient-based technique in 
II20I that selects the achievable rate vector maximizing the sum 
of weighted rates where the weights are marginal utilities at 
each slot, it suffices to solve the following slot-by-slot problem 

'More perfomiance gain may be achieved by canceling inter-cell interfer- 
ence using signal level coordination, such as CoMP (Coordinated Multi Point 
Transmission and Reception) addressed in the LTE-Advanced, which is beyond 
the scope of this paper. 
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(P) which produces the long-term rates that is the optimal 
solution of the (Long-term P). 



(P): max /) = V w,. V r^"(p„ J,) (8) 

p,i ^ — ^ — ^ 

keK. ses 

subject to < 1, WneJV,s eS, (9) 

feG/C„ 

J2pn^pn,max^ Vn S A/", (10) 

VneA/',se5, (11) 



se5 



n ^ pn.mask 



where Wk > is the derivative of its utility Wk = 
'^^Ir'''^ \Rk=Rk(t) that can be interpreted as the weight of user 
k for the slot. For example, we can set Wk as the inverse 
of its average throughput 1/Rk{t) to achieve proportional 
fairness among users ll2Ti . Since the system objective is a 
nonconvex function of the transmit powers and is also tightly 
coupled with integer variables of the scheduling indicators, 
the problem (P) is a mixed-integer nonlinear programming 
(MINLP). Unfortunately, it is known in f22] that even the 
simplified problem, in which user scheduling issue is eliminated 
(i.e., |/C„| = 1 for all n £ A/"), is computationally intractable. 
To find a global optimal solution, we need to fully search the 
space of the feasible powers for all BSs with a small granularity 
along with the all possible combinations of user scheduling. 
Thus, even for a centralized algorithm, it may not be feasible 
in practical systems to solve (P) at each slot. 

III. REFIM: Reference User Based Interference 
Management 

A. Joint Power Allocation and User Scheduling 

User scheduling for fixed power allocation 

We now present our proposed approach to solve the problem 
(P). Note first that for any given feasible power allocation, 
the original problem can be decomposed into intra-cell user 
scheduling problems. 

Lemma 3.1: For any fixed feasible power allocation p, the 
original problem (P) can be reduced to N x S independent 
subproblems for each BS n and subchannel s as follows: 



max 

Is 



subject to 



E 
E 

fce/c„ 



Tk.n 



< 1. 



(12) 
(13) 



Proof: For the given power allocation p, we can rewrite 
h{p, I) as follows: 



Hp, I) 



E E 

n6Af fce/C„ 



Wk^I. 



k.n 



EE 



As Wk and are given parameters, we only have to 

investigate dependencies among Since the constraint (|9]l 
for the given BS n and subchannel s does not affect the 
other BSs and subchannels at all, the original problem can be 



decomposed and is equivalent to individually solving the N xS 
subproblems in (fT2] | and (fT3] i for each BS and subchannel. ■ 

Accordingly, an optimal user scheduling algorithm under the 
given power p is easily obtained by 

rkn f 1: if = Hn, s) = arg max Wk ■ 
1 0, otherwise. 



(14) 



Power allocation for fixed user scheduling 

For any given user scheduling J, the original problem 
reduces to the following power allocation problem: 



max 
p 



l0g2 1- 



k(n,s),n ri 



Vn e A/", s e S. 



subject to ^p^<P"'™ 

n ^ pn^mask 
Ps — s ' 

Solving the above problem requires the knowledge of all 
interference channel gains across cells and noise power, forcing 
it to operate in a centralized fashion. 

To overcome this complexity and develop a distributed 
scheme with small message passing, we introduce the con- 
cept of reference user. Let J\f{n) be the set of neighboring 
BSfl of BS n, and further denote by A{n, s) the set of all 
scheduled users on subchannel s in J\f{n), i.e., A{n, s) ~ 
{k{m, s) I m G J\f{n)}. The reference user for BS n on sub- 
channel s is defined as the user among A{n, s) which has 
the strongest channel gain between the user and BS n on 
subchannel s. We separately denote the index of BS to which 
the reference user belongs and the reference user by ref" 
and fc(ref^,s), respectively. We will elaborate on the way 
of choosing the reference users shortly at the end of this 
subsection. 

Once the reference user is fixed, each BS tries to find its own 
power allocation taking into account just one reference user per 
subchannel instead of solving the above problem considering 
all (N number of) cochannel users in the network. This 
approximation comes from the intuition that just considering 
the worst-case user may be a good approximation of the case 
when all the users are included. The problem (Pn) to be solved 
by each BS n can be written as follows. 



(P„) : max ^ < loga 1- 



ses 



(15) 



log2 1- 



ref'^ ref'^ 

gs "Ps 



ses 

subject to ^p1'<P" 

ses 

p:<p. 



9s " pT + <^s " 



mT^ref" 



n ^ pn.,mask^ ^ g 



(16) 
(17) 



-If there is any chance that users in BS n will handed over to certain 
adjacent BSs, then we consider such a set of BSs as the neighboring BSs of 
BS n, denoted by N{n). This set can be determined a priori by mobile network 
operators at the time of deployment and/or maintained in online based on the 
signal strengths between the BSs. 
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Note that we replace the index of scheduled users k{m,s) 
with the corresponding index of BSs m to keep our notations 
simplell 

For any given user scheduling and reference user selec- 
tion, the corresponding optimal power allocation must satisfy 
Karush-Kuhn-Tucker (KKT) conditions |23]. In particular, let 
A" denote the nonnegative Lagrange multiplier associated with 
the total power budget constraint (flTt . Then, the optimal A" 
and must satisfy the following equalities: 



TABLE I 

General Algorithm Description 



A"ln2 + t" 



ref ^ ref ^ , 



where fl 



I ref"' 
Is 



Eref ,1 I , It 
(GA/- 5^* Ps+ 



pn 



= 0, 



(18) 



(19) 



where [•]" = min [max [•, a] , 6] and jf^" {■) is the received 
SINR as defined earlier in (|4). Note that the term can be 
interpreted as a taxation (or penalty) term. If the reference user 
is close to the BS n (i.e., high interference channel gain glf ° '"), 
then the value of t" increases. Consequently, BS n lowers its 
power level to reduce the harm to the reference user 

Since the modified problem {Pn) is still nonconvex ifTTl , 
lITS I, (fTSl l and fT% are the first-order necessary conditions. 
Therefore, there might exist a duality gap to the optimal 
primal solution. However, encouraged by the state-of-the-art 
asymptotic result [22| that the duality gap becomes zero when 
the number of subchannels is large, we develop an effective 
approximation algorithm for the problem (Pn) based on the 
conditions (fT8]l-(fT9]l. Note that a fixed point equation of in 
(fTSl l is a monotonic function of A". Thus, it can be solved via a 
fast bisection method. Starting from an initial power allocation 
and A", we calculate the power p" and taxation term t" in 
(fTsl l for all subchannels. If the sum of updated power exceeds 
pn,max^ then A" is increased. Otherwise, A" is decreased. With 
the updated power, we repeat this until the equation ([19) holds. 
If no positive value of A" matches the equality, then A" is set 
to be zero. In the latter case (interfering too much), the BS n 
does not use all of its available power. 

Remark 3.2: If the taxation term t" is ignored, our power 
allocation algorithm is reduced to the water-filling (WF) algo- 
rithm ll24l . where each BS acts selfishly in order to maximize 
its own performance. By adding this term > 0, each BS 
operates in a social way by considering the reference user and 
lowers the water-filling level, which could lead to a globally 
better solution. 

General algorithm description: joint user scheduling and 
power allocation 

TABLE U describes a conceptual pseudo-code of the general 
algorithm for our problem. At each slot, each BS starts from a 



'We use the index of BSs instead of the index of scheduled users as 

^ k(rL,s),n 71, m , k(n,s),m „ 

k(n.s) ref"' ref" fcfref"' ,s) ,ref" 

- (7s ' \ Ws ' ^ ^^fc(i-ef^,s), 9s ^ 9s 
ref"" ,m fcfref"" ,s) ,m ref^" kirefV' ,s) ,ref ^ , ref"" fcfref" ,s) 

9s ^ 9s ■ 7s 7s = and = 



follows: 

k(n,s) ,n 
7s ^ ' , cr: 



1: Power initialization 

2: repeat (user scheduling loop): 

3: User scheduling 

4: Neighboring environment abstraction 

5: repeat (power allocation loop): 

6: Update the effect to the neighboring environment 

7: Power allocation 

8: until p converges or max # of iterations is reached 

9: until / converges or max # of iterations is reached 



proper power allocation. For the given power, each BS first ex- 
ecutes the user scheduling and then abstracts what is happening 
in the neighboring environment. For example, the concept of 
reference user can be used as one of the abstraction methods. 
For this given user scheduling and neighboring environment 
abstraction, the BS iteratively updates its power and effects to 
the neighboring environment until p converges or the maximum 
number of iterations is reached. Then each BS repeats the user 
scheduling and the neighboring environment abstraction for the 
updated power and goes into the power allocation loop again. 
This procedure is repeated until the user scheduling / converges 
or the maximum number of iterations is reached. 

The general algorithm not only has a prohibitively high 
computational complexity due to inner and outer loops, but 
also requires multiple information exchanges per slot between 
BSs to reflect the updated interference level followed by the 
updated power. However, the multiple feedbacks in a single slot 
are practically impossible because each user can send its own 
measurement information to the BS once per slot. To overcome 
this complexity, we will propose a simplified algorithm in 
subsection IIII-EI which executes user scheduling and power 
allocation step-by-step without loops. 

B. Online Reference User Selection Method 

Based on the notion of reference user, each BS needs to 
consider the only one user instead of all scheduled users in the 
network on each subchannel. 

From BS n's point of view, although its transmit power 
interferes with all the scheduled users on the corresponding 
subchannel in other BSs, the effect will be critical especially 
for the the user who has the strongest channel gain (or the most 
influenced user from BS n). This dominant victim user can be 
found in the neighboring (or adjacent) BSs J\f{n). Therefore, 
we propose an online reference user selection method, in which 
each BS n chooses the reference user on each subchannel s as 
follows: 

Reference User Selection Rule 

A!(ref",i') IH the reference uyer oti Hubchiitmel 



where ref" 



arg max g 



k{rn,s) .n 



(20) 



It is worthwhile mentioning that the reference user is indepen- 
dently and locally selected by each BS on each subchannel and 
thus no centralized coordination is necessary. Fig. [T]depicts an 
example of our online reference user selection procedure. 
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Step 1. Receive Information about tine Step 3. Solve an approximate per-BS 
scheduled users from neighborign BSs. power allocation optimization problem. 



Fig. 1. Online reference user selection method. 

The rule in ( |20] i provides a guideline for choosing a reference 
user. There may be other methods such as (i) making one 
virtual user by averaging channels of the scheduled users 
from neighboring BSs, and (ii) selecting multiple reference 
users (e.g., select the AI worst users). As we will show later 
in subsection IV-AI such variations do not lead to the high 
performance improvement, compared to the high increase in 
complexity. 

C. Feedback Reduction 

To determine a reference user at each slot, each BS n 
requires (FO) the channel gain g^(™''*)'" of the scheduled 
users in neighboring cells m e Af{n). We call these users 
the candidates of reference users. Once the reference user 
is selected, to calculate a taxation term for power control, 
additional information about the reference user is necessary. 
The followings are the required information for the reference 
user: 

(Fl) The weight of the reference user Wkiref^^s), 

(F2) The received signal strength of the reference user 

fc(ref^ ,s) .ref^^ ref" 

(F3) The noise plus interference strength of the reference user 

The information for the candidates needs to be collected 
by neighboring BSs and be forwarded to BS n. However, 
femto BSs use their Internet connection as backhauls, and 
thus message exchange may not be reliable and fast enough. 
Besides, there is no guarantee on the feedback latency. Even 
for macro BSs with a dedicated backhaul network, the per-slot 
message exchange may be a large overhead. We present a more 
practical solution to reduce the backhaul feedback overhead 
both temporally and spatially. 

Temporal feedback reduction: Instead of the per-slot message 
exchange for all the information, (i) each candidate user first 
calculates the time-average of the information and send them 
to its associated BS infrequently (say, every T ^ 1 slots), (ii) 
and then the BS broadcasts the information about all candidate 
users to its neighboring BSs through wired backhauls. The only 
thing that the BSs exchange at each slot is the indexes of the 
scheduled users. Each BS will maintain a table that contains 
these averages of candidate users. Once each BS receives the 
indexes of scheduled users from neighboring BSs, then it uses 



the information in the table corresponding to the indexes. Note 
that the slow feedback mechanism can be applied to femto BSs. 

Spatial feedback reduction: First, we reduce the amount of 
infrequent feedback by making macro BSs send the information 
only for edge users. This idea comes from the intuition that the 
users in the center of cells are not likely to be selected as the 
reference user because the center users do not receive too much 
interference. 

Second, we eliminate the per-slot feedback of indexes of 
scheduled users for femto BSs. The indexes of scheduled users 
(very small amount of information) can be easily exchanged 
between neighboring macro BSs at each slot through dedicated 
backhauls. However, this per-slot message exchange is not pos- 
sible for femto BSs because their backhauls do not provide any 
guarantee on the feedback latency. To overcome this difficulty, 
we propose an alternative solution for the femto BSs that does 
not require per-slot message exchange at all. In the proposed 
solution, the femto BSs obtain the indexes of scheduled users 
by overhearing downlink control message (e.g, DL-map in the 
IEEE 802. 16e |25|) from neighboring macro BS. This may 
need slight modification frame structure for the femto BS in 
the current wireless standards. 

The remaining challenge is on the reverse direction, i.e., 
sending the indexes of scheduled users in the femto BSs to 
the neighboring macro BSs. We pay attention to one of the 
main features of femto cells, that is, the small coverage. In 
other words, the distances between users in a femto BS are 
relatively much shorter than the distances from the neighboring 
macro BSs. Thus, from neighboring macro BSs' perspective, it 
seems reasonable to assume that the users in the femto cell 
are spatially located at the same point. Based on this spatial 
simplification, the neighboring macro BSs simply can pick any 
user in the femto cell and consider as the candidate user. 

D. Initial Power Setting 

Our algorithm requires an initial power value to compute 
the power allocation (see line 1 in Table I). Since our problem 
is a nonconvex problem, different starting points may lead 
to different solutions with different speeds. The following 
three strategies for the choice of initial power are carefully 
investigated in this paper 

1) Uniform rule. The power allocation always starts from 
the same point for every slot. Each BS uniformly spUts 
its maximum transmission power to all subchannels, i.e., 

^n^init^'^ pn.max j g 

2) Random rule. Each BS randomly chooses the initial power 
level for each subchannel between and pn.max ^j^^ 
then each BS scales it up with an appropriate weight 

to use up the total transmission 
power budget, i.e., Es-P"'™*(0 = P"'™"^. 

3) Previous rule. Each BS starts from the power used at the 
previous slot, i.e., — p'^[t — 1). 

We will demonstrate later in subsection IV-AI that although we 
avoids multiple loops in a slot for power allocation and user 
scheduling, the performance of the previous rule is likely to 
remain unchanged while the other rules lose the performance 
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TABLE II 

Complexity comparison of various algorithms 



Algorithms 


Computational complexity 


Signaling complexity (inter-BSs) 




User scheduling 


Power allocation 


Per-slot feedback 


Periodic feedback 


EQ 


0(SK) 


Zero 


Zero 


Zero 


WF 


0{SK) 




Zero 


Zero 


MGR 


0(SK) 


-^0{S)+n„0{SK) 


Zero 


\N^\S 


REFIM 


0(SK) 




macro: pS & femto: Zero 


p\K^\AS 


MC-IIWF 


Ti ■ {0{SK) + 0{SN)) 


Complete infoiTnation is assumed. 



much. This is because in some sense the previous rule exploit- 
ing temporal correlation has the effect of iterations for power 
allocation in a slot-by-slot manner. This result encourages us 
to design a simplified algorithm in subsection IIII-EI which gets 
rid of multiple loops in a slot and executes user scheduling and 
power allocation sequentially with the previous rule. 

E. REFIM: Reference Based Interference Management 

We now propose our final algorithm, called REFIM (REFer- 
ence based Interference Management), in Table |III] that merges 
all the components developed in the above. REFIM adopts the 
previous rule for initial power setting (see line 1), and uses 
the notion of reference user for the neighboring environment 
abstraction and limits the number of reference users to one ( see 
line 3). While the general algorithm has user scheduling and 
power allocation loops (see lines 2 and 5 in Table |l]i, REFIM 
executes user scheduling (see line 2) and power allocation (see 
line 5) sequentially without loops. This step-by-step approach 
can not only be done very fast in a slot, but it requires the 
feedback from each user just once per slot. Through extensive 
simulations, such a simple algorithm will be shown to be 
efficient. 

TABLE III 

REFIM: Reference based Interference Management 



1: Power initialization p"(t) p"(t — 1) 
2: User scheduling according to ([14): 

k{n,s) = arg max uij. • r^'"{p"). 

3: Reference user selection according to 1 1201 . 
4: Taxation update according to l |18t : 
5: Power allocation via bisection: 
[a,b] 4- [0,A^,J. 
while \ J2sPs - P"'"""^! < <5, 

Set A" = (a + b)/2 and update p" according to )18t . 
if T,sPs > P"'""^, then [a,b] <- [A",fe], 
else f^sPs < P"'™"', then [a,b] ^ [a, A"], 
end while 



IV. Complexity Analysis 

In this section, we analyze the computational complexity 
and inter-BSs signaling complexity of REFIM compared to 
conventional equal power allocation (EQ) and selfish water- 
filling (WF), as well as MGR (Multi-sector GRadient) [I2| and 
MC-IIWF (MultiCell Improved Iterative Water Filling) ifTTl . 
Table |ll] summarizes the results of complexity analysis. 



Computational complexity consists of two parts: the com- 
plexity from user scheduling and power allocation. User 
scheduling has a linear complexity 0{SK) with the number 
of users for each subchannel for all algorithms except MC- 
IIWF. For power allocation, EQ has zero complexity. WF 
can be obtained by either an exact algorithm 0{S) or an 
iterative algorithm (i.e., a bisection method that converges to a 
solution with a certain error tolerance e) O (log2 ^"1"'^ ) Ii24l . 
The only difference between WF and REFIM is the taxation 
term considering the reference user Thus, the complexity of 
power allocation for REFIM is basically the same as that for 
WF. MGR in |[l2l, one of the state-of-the-art dynamic IM 
algorithms, adjusts the power allocation for every rip > 1 slots 
(relatively slowly) and condenses the complexity for updating 
power to l/up. However, MGR has additional complexity from 
a virtual scheduling that needs to be run riy times per slot. Ac- 
cordingly, the total computational complexity for power alloca- 
tion is high, -^0{S) + nyO{SK). Another recently developed 
MC-IIWF in 1 11) is the centralized algorithm that has iteration 
loops for power allocation and user scheduling. Let Ti be the 
number of iterations needed for iteration loops. Then, the total 
computational complexity is equal to Ti • {0{SK) + 0{SN)). 

Now let us investigate the inter-BSs signaling complexity. EQ 
and WF do not require any inter-BS message passing overhead 
because they are autonomous algorithms without considering 
neighboring BSs, but at the cost of performance reduction, as 
will be shown in Section |Vl MGR adjusts the power allocation 
slowly so that it requires not per-slot but periodic feedback, 
jA^nlS* (sensitivity information for neighboring BSs and all 
subchannels). MC-IIWF, a centralized algorithm, assumes a 
central control unit to have complete information. REFIM 
requires the periodic feedback about the candidate users for the 
reference users, p|/C„|A5, where p is the average percentage 
of edge users and A = 4 is the number of required information 
(F0)~(F3) about the reference user Note that while macro BSs 
require the per-slot feedback for the indexes of scheduled users 
at each subchannel, femto BSs do not. 

In summary, the computational complexity of REFIM is the 
same as that of WF and is much lower than that of state-of- 
the-art dynamic IM algorithms such as MGR and MC-IIWF. 
For signaling complexity, although the feedback per slot-wise 
manner is challenging, the information needs to be exchanged 
only between neighboring macro BSs at each slot are just the 
indexes of scheduled users. We believe that such small amount 
of information can be easily exchanged through high-speed 
dedicated backbones. 
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V. Performance Evaluation 

We verify the system performance through extensive simula- 
tions under various topologies and scenarios. First, in order to 
verify the effectiveness of REFIM by varying several tunable 
parameters and comparing it with other algorithms, a two- 
tier macro-cell networks composed of 19 hexagonal cells is 
considered in subsection I V-AI Second, in order to provide more 
realistic simulation results, a real 3G BS deployment topology 
consisting of heterogeneous environments (urban, suburban 
and rural areas) is considered in subsection IV-CI Third, a 
heterogeneous network topology with small cells inside macro 
cells is also considered in subsection IV-DI 

In our simulations, macro and small cells are loaded with 20 
and 4 users, respectively and they are uniformly distributed in 
each cell. All users are assumed to have a logarithmic utility 
function, i.e., U{Rk) — log Rk, but the other utility functions 
enforcing more fairness are also considered in our technical 
report [26|. We consider a system having 16 subchannels 
each of which consists of multiple subcarriers. The maximum 
transmit powers of macro and small BSs are 43dBm and 15dBm 
127 1. respectively. In modeling the propagation environment, an 
ITU PED-B path loss model 16.62+37.6 log]^Q((i[m]) for macro 
cells and an indoor path loss model 37 + 32\ogiQ{d[Tn]) for 
small cells with lOdB penetration loss due to walls are adopted. 
Jakes' Rayleigh model (with the speed of 3km/h), where the 
channel coefficients vary continuously over slots, is adopted for 
fast fading. The channel bandwidth and the time-slot length are 
set to be lOMHz and 1ms, respectively. 

REFIM is compared to conventional EQ and WF, as well 
as MGR L12J and MC-IIWF [11] developed recently. As 
performance metrics, the geometric average of user throughputs 
(GAT) and the average of edge user throughputs (AET) are 
used. We use GAT since maximizing this metric is equivalent 
to our system objective (sum of log throughputs). We consider 
AET as the average of the bottom 5% of user throughputs 
(i.e., 5th percentile throughput), which can be regarded as a 
representative performance metric of cell-edge users. 

A. Ejfectiveness of Proposed Algorithm 

Fig. |2(a)| shows the GAT performance of REFIM for the 
number of reference users per subchannel and that of EQ 
as a baseline. As mentioned in Remark 13.21 REFIM without 
a reference user is reduced as WF. REFIM taking reference 
users into consideration can obtain higher performance gain. 
It is noteworthy that considering only the one reference user 
per subchannel is efficient enough because it can obtain more 
than 97% of the performance considering all the six neighbors. 
Fig. |2(b)| shows the effect of iteration loops and initial power: 
(i) adding user scheduling loop and/or power allocation loop 
gives additional performance gain from any initial power setting 
and (ii) using power at the previous slot as an initial power 
outperforms other two strategies. However, for the case in 
which power at the previous slot is used as an initial power, 
the performance gain from adding power allocation loop is 
marginal. Fig. |2(b)| is an encouraging result that leads us to 
design the algorithm without loops and use the previous power 
as an initial power 
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<(C— WF 
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Number of reference users per subchannel 

(a) The number of reference users 




(X,X) (X,0) (0,0) 

(user scheduling loop, power allocation loop) 

(b) Iteration loops and initial power 



Fig. 2. Effect of several tunable parameters. 



We conjecture that this is because in some sense using 
the previous power rule has the effect of iterations not in 
a slot but over a series of multiple slots. To verify this 
statement, we provide additional simulations in a Unear two- 
cell network where the distance between BSs is 2km. Each 
cell has two groups of users: center and edge users, whose 
distances from the associated BS are within 200-400m and 700- 
900m, respectively. We adopt the simple network configuration 
to gain insight easily as well as for easy of presentation, but 
all discussions can be extended to the general cases. 

Fig. |3(a)| shows the time-series of transmit powers on dif- 
ferent subchannels. Although the powers do not seem to quite 
converge, they remain the certain levels for several dozens of 
slots. In Fig. |3(b)[ we plot the average transmit power levels 
on different subchannels during the period between 75 and 125 
slots. As can be clearly seen, each BS exclusively utilizes five 
subchannels of high powers and shares six subchannels of low 
powers with the other We further investigate the relationship 
between the user groups and their subchannels on which they 
are scheduled. Interestingly, for the most of time (more than 
98% of slots), the center and edge users are served by the set 
of subchannels with low and high powers, respectively. In other 
words, it is highly likely that on each subchannel a user will be 
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(a) Time series of transmit powers for different subchannels 

4r 
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a 2.5 




2 
1.5 
1 
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Subchannel index 

(b) Average transmit powers between 75 and 125 slots 

Fig. 3. Transmit power levels for different subchannels in a Hnear two-cell 
network. 



selected by the scheduler who has similar channel conditions to 
the user selected at the previous slot. Thus, given similar user 
scheduling over consecutive slots, the power allocation with 
the previous power can be interpreted as if it has iterations 
over a series of multiple slots. This explains why REFIM that 
executes user scheduling and power allocation step-by-step in 
a slot can achieve a good solution without much performance 
degradation. 

In Fig. m we also test the effect of outdated feedback 
information about reference users. We consider two types 
of users with different speeds: nomadic users (or stationary 
users that have fixed path loss and shadowing factor, but have 
time-varying fast-fading) and mobile users (moving fast with 
speed of 60km/h). For nomadic users, the GAT performance 
degradation is relatively small, even though we choose a 
long feedback period such as 200 slots. For mobile users, 
the GAT performance naturally decreases due to the error 
of feedback information as the feedback period increases. In 
practical systems, different types of users with different speeds 
are expected to coexist. In such an environment, to reduce 
the amount of feedbacks while maintaining the performance 
degradation marginal, it is essential to adaptively control the 
period of feedback for the different speed users, e.g., T — 200 




20 40 60 80 100 120 140 160 180 200 
Period of feedbacl< 

Fig. 4. The effect of outdated feedback information. 
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Fig. 5. Comparison with other algorithms. 



slots for nomadic users ('■ 
users (-60km/hfl 



^3km/h) and T = 10 slots for mobile 



B. Performance Comparison with Other Algorithms 

Now we compare the performance of REFIM with other four 
algorithms: conventional EQ, selfish WF, MGR [12| which 
adjusts power allocation infrequently, compared to per-slot 
basis user scheduling, and MC-IIWF 1 1 1 1 which is a centralized 
algorithm achieving the near-optimal performance. Fig. |5]shows 
the CDF (cumulative distribution function) of the throughput of 
entire users in the network for different algorithms. Compared 
to EQ, WF and MGR, REFIM can improve the throughputs 
for all users in the network. Particularly, we can observe 
higher improvement (43% improvement in AET compared to 
EQ) for users achieving low throughputs, i.e., users at cell 
edges. This is due to the fact that IM is mainly targeted 
for performance improvement of cell-edge users. In addition, 
REFIM can achieve about 95% of the performance of near- 
optimal MC-IIWF in terms of two representative throughput 
metrics (GAT and AET) as well as the arithmetic average of 

''We use the simplified version of random waypoint model (28 1, where 
each user starts from an initial point, randomly chooses its destination, and 
moves toward it at a given speed. After reaching the destination, it repeats this 
process unite the end of simulation time. 
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Fig. 6. Real 3G BS deployment map. 



user throughputs (AAT). It is somewhat surprising that such a 
simple distributed algorithm can obtain a similar performance 
to the centralized algorithm that is hard to implement due to 
prohibitive complexity. 

C. Topology with Real BS Deployment: Urban, Suburban and 
Rural Environments 

Fig. |6] depicts the map of BS layout that we use for more 
realistic simulations. It is a part of real 3G network operated 
by one of the major mobile network operators in Korea. There 
are a total of 30 BSs within 20 x 10 km^ rectangular area. We 
assume that the number of BSs per unit area is proportional to 
the user density. In other words, the average number of users 
per cell is almost similar because BSs in an urban environment 
cover a small area and BSs in a rural environment a large 
area. Under this assumption, we generate users one-by-one 
in the rectangular area and attach them to the closest BS 
until each BS will have 20 users. We choose this partial map 
to include a challenging scenario that three environments are 
mixed together We tested several other maps, and obtained 
similar or even better performance of REFIM. 

We examine three different zones: urban (15 BSs in 4.5 x 
4.5 km^), suburban (15 BSs in 12 x 6 km^) and rural (8 BSs 
in 9x9 km^) area^ Fig. |7] shows GAT and AET performance 
under urban, suburban and rural environments. As expected, 
we can obtain high performance improvements in the urban 
and suburban environments. However, almost low or no gain 
is found in the rural environment, which means that the IM 
does not take much effect in a sparse topology, which follows 
our intuition. For example, in the urban area, the performance 
gains of REFIM are 16% and 42% in terms of GAT and AET, 
respectively. 

Another nice feature of REFIM is incremental deployment. 
Suppose that we implement our algorithm only on the BSs in 
a specific area. While the BSs inside this area performs well as 
we want, the BSs in the boundary of the area does not. This is 
because they may not receive information about reference users 
from the some of its neighboring BSs on which our algorithm 
is not implemented. Even in such a case, our algorithm will 

^In order to see clearly how the density of BSs affects the performance 
gain, the distance between BSs in suburban and rural zones are increased by 
1.5 and 2 times, respectively. 
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Fig. 7. Throughput performances in real BS deployment topology. 
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Fig. 8. Effect of partial deployment. 



automatically reduce to WF. Thus, it performs like WF at least 
and better than EQ. Compared to the full deployment case, the 
partial deployment case where only 15 BSs (mainly selected 
from the urban areas among 30 BSs) are equipped with REFIM 
can achieve more than 85% gain as shown in Fig. [8] The result 
encourages the mobile network operator to upgrade its BSs 
incrementally from the urgent ones, e.g., densely located BSs 
experiencing heavy interference, and thus low capacity. 
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Fig. 9. Performance in the heterogeneous network topology. 



D. Heterogeneous Networks (Macro + Small BSs) 

Now we consider a heterogeneous network topology having 
several small cells inside macro cells as shown in Fig. |9(a)| 
As a kind of small BSs, we consider the femto BSs that are 
deployed and provide a high-speed indoor access mainly to 
home users. We need to reflect femto-to-femto cell interference 
as well as macro-to-macro and macro-to-femto cell interfer- 
ences. To this end, we consider the mixture of three femto BSs 
deployment cases: (i) only one femto BS case (no femto-to- 
femto interference), (ii) two symmetric femto BSs case (strong 
femto-to-femto interference): two femto cells are adjacent with 
each other and each femto BS is located in the center of the 
home, (iii) two asymmetric femto BSs case (very strong femto- 
to-femto interference): two femto cells are adjacent with each 
other and femto BS 1 in home 1 is located at the border between 
homes. The last case can often happen because users locate their 
own femto BSs wherever they want without considering next 
door neighbors. 

In Figs. |9(b)| and |9(c)[ we compare the performance of 
the spectrum sharing policy (i.e., universal frequency reuse) 
between macro and femto cells with that of the spectrum 
splitting policy where macro and femto cells orthogonally use 



the resource. Note that the performance curves for the spectrum 
sharing policy and the spectrum splitting policy are represented 
by solid and dotted lines, respectively. The x-axis represents 
the ratio of subchannels used by macro cells among all 16 
subchannels. 

In the case of five femto cells per a macro cell in Fig. 
|9(b)| there exists small cross-tier (macro-to-femto) interference. 
Thus, the spectrum sharing policy is always better than the 
spectrum splitting policy. For example, even the performance 
of HQ without any interference management in the spectrum 
sharing policy is higher than or equal to the performance 
of REFIM in the optimal spectrum splitting (at 8/16). If the 
number of femto cells increases, then the portion of macro 
users who will see more and closer femto cells increases. 
Consequently, it is highly probable that their performances are 
degraded by the severe cross-tier interference. As shown in Fig. 
|9(c)| with ten femto cells per a macro cell, the best performance 
in the optimal spectrum splitting (at 6/16) can catch up with the 
that of HQ and WF in the spectrum sharing policy. However, 
if the proposed IM algorithm, e.g., REFIM, is adopted in the 
spectrum sharing policy, then we can mitigate the cross-tier 
interference, resulting in the better performance than any case 
in the spectrum splitting policy. Note that the performance 
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improvements of our REFIM in the spectrum sharing poUcy 
compared to EQ in the spectrum sharing and spectrum splitting 
policies are 38% and 106% in Fig. |9(b)| and they become larger 
as the number of femto cells increases, i.e., 56% and 107% in 
Fig.|9(^ 

The spectrum splitting policy is expected to be widely used 
rather than the spectrum sharing poUcy in an early stage of 
femto cell deployments because it makes the femto cells easily 
coexist with macro cells without worrying about the macro-to 
femto interference. However, the results in Fig. |9] enlighten us 
on the potential gain of the spectrum sharing policy. Therefore, 
we believe that, if the IM algorithms become more mature in the 
near future, then the spectrum sharing policy will be adopted 
in order to maximally exploit the spectral resources due to the 
explosive traffic demands 

VI. Conclusion 

Heterogeneous access networks, consisting of cells with 
different sizes and ranging from macro to femto cells, will 
play a pivotal role in the next-generation broadband wireless 
network. They can increase the network capacity significantly 
to meet the explosive traffic demand of users with limited 
capital/operating expenditures and spectrum constraints. One 
of the biggest challenges in such environments is how to 
effectively manage interferences between heterogeneous cells. 
To tackle this challenge, this paper developed REFIM, which is 
an efficient low-complex and fully distributed IM in downlink 
heterogeneous multi-cell networks. Our key idea is to use the 
notion of reference user, which can spatially simplify the impact 
of all other neighboring cells by a single virtual user and result 
in the power control algorithm with low computational and 
signaling overhead. In order for REFIM to be implemented even 
on the femto BSs, we also further reduced the feedback over 
backhauls both temporally and spatially. Through extensive 
simulations and complexity analysis, we demonstrated that 
REFIM not only performs well but also is practically imple- 
mentable. We also concluded that as long as appropriate IM 
algorithms such as REFIM are adopted, the spectrum sharing 
policy can outperform the best spectrum splitting policy where 
the number of subchannels is optimally divided between macro 
and femto cells. 
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